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ABSTRACT 

Aims. We investigate the temporal evolution of magnetic flux emerging within a granule in the quiet-Sun internetwork at disk center. 
Methods. We combined IR spectropolarimetry of high angular resolution performed in two Fei lines at 1565 nm with speckle- 
reconstructed G-band imaging. We determined the magnetic field parameters by a LTE inversion of the full Stokes vector using the 
SIR code, and followed their evolution in time. To interpret the observations, we created a geometrical model of a rising loop in 3D. 
The relevant parameters of the loop were matched to the observations where possible. We then synthesized spectra from the 3D model 
for a comparison to the observations. 

Results. We found signatures of magnetic flux emergence within a growing granule. In the early phases, a horizontal magnetic field 
with a distinct linear polarization signal dominated the emerging flux. Later on, two patches of opposite circular polarization signal 
appeared symmetrically on either side of the linear polarization patch, indicating a small loop-like structure. The mean magnetic flux 
density of this loop was roughly 450 G, with a total magnetic flux of around 3xl0 17 Mx. During the ~12min episode of loop occur- 
rence, the spatial extent of the loop increased from about 1 to 2 arcsec. The middle part of the appearing feature was blueshifted during 
its occurrence, supporting the scenario of an emerging loop. There is also clear evidence for the interaction of one loop footpoint with 
a preexisting magnetic structure of opposite polarity. The temporal evolution of the observed spectra is reproduced to first order by 
the spectra derived from the geometrical model. During the phase of clearest visibility of the loop in the observations, the observed 
and synthetic spectra match quantitatively. 

Conclusions. The observed event can be explained as a case of flux emergence in the shape of a small-scale loop. The fast disappear- 
ance of the loop at the end could possibly be due to magnetic reconnection. 
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1. Introduction 

Our knowledge of small-scale magnetic fields located in the 
photosphere and low chromosphere has evolved extremely fast 
thanks to the high spatial resolution observations done with the 
Hinode satellite (e.g., Lites et al. 2008), and the improvement 
of the spatial resolution of large ground-based solar telescopes 
achieved by adaptive optics systems (e.g., von der Liihe et al. 
2003 ). The results of such observations tracing the weakest mag- 
netic fields can now be compared directly with theoretical mod- 
els, for example on the question whether a local dynamo oper- 
ates in the solar photosphere (Cattaneo 1999 ). The observations 
may also help to refine modern numerical simulations of sur- 
face magnetoconvection (Vogler et al. 120051 Schaffenberger et 
al. 120061 Stein & Nordlund [20061 Abbett [20071 by providing 
more stringent boundary conditions on magnetic field properties 
in the solar photosphere. 

Direct measurements of the 3D topology and the tempo- 
ral evolution of magnetic loops associated with small granular 
and intergranular structures are, however, still rare. The mea- 
surements require accurate spectropolarimetric 2D data of all 
Stokes parameters with high spatial resolution, high signal-to- 
noise ratio, and high cadence. Martinez Gonzalez et al. (J2007 ) 
reported observations of the magnetic field vector in the internet- 
work. They showed evidence that small-scale (only 2-6 arcsec 



long) low-lying loops connect at least 10-12% of the internet- 
work magnetic flux measured at 1565 nm. 

Centeno et al. (120071 CE07) and Orozco Suarez et al. d2008l 
OR08) analyzed the temporal evolution of such small-scale mag- 
netic features. They used similar data sets of the Fei630.15nm 
and 630.25 nm spectral lines from the Hinode/SOT spectropo- 
larimeter (Lites et al. I200T1 Kosugi et al. |2"U07l ). OR08 con- 
structed maps of the temporal evolution of linear and circular 
polarization signals, line-of-sight velocity, and the intensities of 
the Fe 1 630.25 nm and Ca n H spectral line. They found magnetic 
signals appearing at the central parts of granules, but without a 
significant linear polarization signal. Estimated lifetimes of the 
events were 20 and 14min. They did not find any significant 
coupling of the events to the chromosphere. CE07 applied an in- 
version assuming local thermodynamic equilibrium (LTE) to the 
Hinode spectra to retrieve information on the temporal evolution 
of the magnetic flux and its topology. Similar to OR08, they also 
reported the appearance of magnetic signals in the central parts 
of granules, and documented a drift of vertical dipoles toward the 
surrounding intergranular lanes. CE07 found a slightly shorter 
lifetime of only 8 min. Moreover, they found that the appearance 
of horizontal magnetic fields precedes that of the vertical fields. 
They interpreted their results as the rise of small-scale magnetic 
loops through the photosphere, although they could not exclude 
a descending loop. 
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In this paper, we present a similar analysis of the tempo- 
ral evolution of the appearance and topology changes during a 
particular flux emergence event in the internetwork at the disk 
center, using ground-based observations in the infrared spectral 
region rendering high polarimetric precision and magnetic sen- 
sitivity. We compare these observations with synthetic spectra 
from a geometrical model of a rising loop to investigate whether 
the observations can be interpreted in such terms. 

2. Observations 

We observed a very quiet region at disk center using the German 
Vacuum Tower Telescope (VTT, Schroter et al. 119851 1 located 
at the Observatorio del Teide (Tenerife, Spain) on June 6, 2008 
between 09:57:50 - 1 1 : 10:30 UT. 

We used the Kiepenheuer-Institut Adaptive Optics System 
(KAOS, Berkefeld 2006) to improve the spatial resolution by a 
real-time correction of wavefront distortions. The uncorrected 
Fried parameter ro reached an average value of about 9 cm dur- 
ing the whole observing run that is sufficient to reach the diffrac- 
tion limit of the VTT of around 076 at infrared (IR) wavelengths 
after the correction. 

The light coming from the telescope was separated into vis- 
ible and IR radiation components using a dichroic beam split- 
ter (splitting edge at 800 nm). The IR radiation was sent to the 
Tenerife Infrared Polarimeter (TIP 2, Collados et al. 2007]) at- 
tached to the main spectrograph of the VTT. The visible frac- 
tion of the radiation was reflected to an optical bench where 
the following instruments were fed in parallel by two additional 
achromatic beam splitters: a G-band broad-band filter (center 
and FWHM: 430.6 nm, 0.9 nm), a narrow-band Zeiss Lyot filter 
for Ha (FWHM -0.025 nm), and the TESOS 2D spectrometer 
(Kentischer et al. [T998 Tritschler et al. |2002| l. More informa- 
tion on the instrumental setup can be found also in Kucera et al. 
( 120081 or Beck et al. (120071 ). 

We used the TIP spectropolarimeter with a 0745 wide and 
77" long slit to repeatedly scan a narrow 5" wide area in 10 
steps. TIP covered the spectral region around the Fe i lines at 
1564.85 nm and 1565.28 nm. For a more detailed description of 
these lines see, e.g., Bellot Rubio et al. (120041) or Khomenko & 
Collados (2007). Spatial and spectral sampling of the data were 
0718/px and 14.64 mA/px, respectively. We acquired the four 
modulation states employed by TIP with 250 ms exposures each. 
To obtain a higher signal-to-noise ratio in the recorded Stokes 
parameters (/, Q, U, V), we then added 8 sets of the modula- 
tion states at each scan position, yielding a total integration time 
of 8 s and a cadence of ~10 s for two consecutive scanning posi- 
tions. The acquisition of the whole 10-step raster thus took about 
100 s for every repetition of the scanning. 

The G-band images were acquired with a 14-bit PCO 4000 
digital camera. We set the exposure time to 5 ms and the cadence 
to 0.6 s; the field of view (FoV) was 70" x 100", and the spatial 
sampling was 0708 per pixel. The Ha data (see, e.g., Kucera et 
al. [2008 ) and the data taken with the TESOS instrument were 
not used for this publication. 

3. Data reduction and analysis 

We photometrically corrected all data by subtracting the dark 
current and applying flat-field frames. To compensate for the 
instrumental polarization, we applied the correction described 
by Schlichenmaier & Collados (2002) to the spectra of TIP, us- 
ing a data set obtained with the IR instrument calibration unit 




I: 0.96-1.04; v: -1.1-1.1 [km/s]; L: 12x10 4 -14.5x10- 4 ; |V|: 9x1 0'-l 5.5x1 0' 4 



Fig. 1. Overview of the TIP dataset. The panels show from top to 
bottom 2D maps of the continuum intensity, the Doppler shifts of 
the Fe 1 1564.85 nm line (black to mid-gray colors correspond to 
the blueshifts, mid-gray to white colors represent the redshifts), 
the net linear polarization, the circular polarization, and the po- 
larity of the detected magnetic structures. The size of a single 
map here is 5"x20"; their corresponding positions inside the G- 
band FoV are marked in Fig. [3] 



of the VTT. Then, for a complete correction, the coelostat con- 
figuration of the telescope was modeled according to Collados 
( 119991 ) and Beck et al. (120051) . We finally estima ted statistically 
the residual crosstalk following Collados (2003). An overview 
of the relevant section of the dataset obtained with TIP is shown 
in Fig.Q] It shows from top to bottom the continuum intensity I c , 
the line-core Doppler shifts of the Fe i 1564.8 nm line, the total 
linear polarization L — J (Q 2 + U 2 ) 1 ^ 2 1 I c dA, the total circular 

polarization | V\ — J \V\/I c dA, and the magnetic field polarity de- 
rived from the order of minimum and maximum of the V signal. 
The region of the flux emergence is marked with white rectan- 
gles; the repetitions of the scanning marked with the consecutive 
numbers 1 to 6 are used for comparison with the geometrical 
loop model later on (see also Fig.|4|. 

To derive the magnetic field parameters, we used the SIR 
code (Stokes Inversion based on Response functions, Ruiz Cobo 
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Fig. 2. Examples of the typical Stokes profiles detected in the middle part (top row) and southern footpoint (bottom row) of the 
emerging magnetic feature. Thin lines show observed profiles, thick lines the best-fit profiles of the inversion. Narrow panels show 
their residual differences. The x, y, t parameters represent spatial positions of extracted profiles and the times of particular exposures. 



& del Toro Iniesta 1992) for a simultaneous inversion of the 
Stokes profiles of Fei lines at 1565 nm in a similar way as de- 
scribed by Bellot Rubio et al. (2004). We used a two-component 
model of the solar atmosphere with a magnetic and a field-free 
component that allows us to derive accurate parameters also in 
the case of spatially unresolved magnetic fields. All magnetic 
parameters were kept constant with height during the inversion. 
To demonstrate the quality of the performed inversion, we dis- 
play the typical Stokes profiles detected in the middle part and in 
the southern footpoint of the emerging feature in Fig.|2]together 
with the corresponding best-fit profiles determined by SIR, and 
their residual differences. The rms noise in continuum windows 
of the observed Stokes profiles shown in Fig.|2]was 3.5 x 10~ 3 , 
1.1 x 10~ 3 , 1.1 x 10~ 3 , and 1 x 10~ 3 for the middle part and 
3.7 x 10~ 3 , 1.1 x 10~ 3 , 1.5 x 10~ 3 , and 8 x 10~ 4 for the southern 
footpoint in the 7, Q, U, and V Stokes profiles, respectively. 

We used the Kiepenheuer-Institut Speckle Interferometry 
Package (KISIP) code (Woger et al. 2008 ) for the speckle re- 
construction of the G-band data. A burst of 100 images was se- 
lected for each reconstruction leading to a temporal cadence of 
63 s. The maximal time difference between a TIP scan position 
and the G-band reconstruction closest in time is thus around 30 
seconds. A sample of the acquired images is given in Fig. [3] 

Because of the importance of a precise co-alignment of all 
data, we took additional grid target data with all involved in- 
struments. Using these data, we determined the rotation between 
the TIP spectropolarimetric rasters and the G-band images. The 
overall shift of these datasets, which also includes the effects of 
differential refraction, was then derived with sub-arcsec preci- 
sion using a cross-correlation technique. At the end, we com- 
pared the co-spatiality of some G-band bright points with en- 
hanced Stokes V signals as an independent check of the applied 
co-alignment method. 



4. Emerging loop in observations and inversion 

The temporal evolution of the polarization signals in the small 
photospheric internetwork area under study, and the magnetic 
field parameters (i.e., magnetic flux density, field inclination, 
and field azimuth) derived by the inversion are shown magni- 



fied in Fig. [4] The first row displays the speckle-reconstructed 
G-band images. The superimposed contours mark areas where 
a significant polarization signal above the noise level was de- 
tected. The green contours mark areas of significant net lin- 
ear polarization. The red and blue contours represent negative 
and positive circular polarization, respectively, given as CP = 
J(|V| X polarity) I I c d A. 

In the first snapshot taken at 10:37:06 UT, three areas of pos- 
itive circular polarization are visible. These elements represent 
the preexisting magnetic field seen in the data set well before, 
during, and after the time period investigated here. The second 
snapshot (10:38:47 UT) shows a new, predominantly linearly po- 
larized signal at (x,y)=(-2", 2"), emerging at the location of a 
granule. In the next minutes, the topology of this feature changes 
significantly. It expands rapidly, but mostly within the granu- 
lar boundaries. The next snapshots taken at 10:42: 10UT and 
10:43:53 UT already show areas of opposite circular polariza- 
tion signal occurring symmetrically on both sides of the region 
with enhanced linear polarization signal. The linear and circu- 
lar polarization patches partly overlap at these moments. This 
configuration could correspond to a small magnetic loop-like 
structure whose footpoints are co-spatial to the areas of opposite 
circular polarization. The distance of the potential footpoints is 
roughly 1-1.5". A slight downward deformation of the preex- 
isting positive circular polarization patch at (-2", 4"), clearly 
already visible at 10:43:53 UT, suggests the onset of an interac- 
tion of the emerging negative polarization patch with this fea- 
ture. Later on, at 10:45:33 UT, the sudden disappearance of the 
linear polarization signal is followed by a weakening of both 
circular polarization patches which have a distance of about 
2" at this moment. Finally, at 10:47: 15 UT and subsequently at 
10:48:57 UT, the remnants of the circularly polarized features 
have disappeared completely, and there is no further evidence of 
their re-occurrence later on in the data set. The whole lifetime of 
the detected magnetic feature thus was less than 12 minutes. 

The overall course of the event seems to be rather asym- 
metric in time. Figure|4] indicates a gradual growth of the area 
of enhanced net linear polarization and the polarization ampli- 
tude itself in the four-panel sequence between 10:38:47 UT and 
10:43:53 UT, followed by a sudden disappearance in the next 
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Fig. 3. Sample of the speckle-reconstructed G-band images showing the surroundings of the observed region. The locations of the 
TIP scans shown in Fig.Q]are marked by white rectangles. The zero point of the coordinate system is shifted 55 arcsec from the disk 
center. 



panel. The flux emergence leaves no detectable brightness pat- 
terns on the underlying granule, which grows steadily in size. 
The emergent magnetic flux also does not influence the appear- 
ance and brightness of G-band bright points in the vicinity of the 
granule significantly; most of them stay visible at about the same 
location for the whole time span of the flux emergence. 

The second row of Fig. [4] shows the temporal evolution of 
the magnetic flux density. Compared with the preexisting ambi- 
ent magnetism of kG strength, the flux density of the emerging 
magnetic feature is much weaker and reaches a mean value of 
about 450 G throughout its 12-min period of occurrence. The 
temporal evolution of this parameter on spatial cuts connecting 
the two footpoints is also shown in the top panel of Fig. [5] The 
displayed error bars correspond to the error estimate of the SIR 
code. Moreover, using the inversion results on the magnetic flux 
density and the magnetic filling fraction, we estimated the total 
magnetic flux seen in each of the two footpoints to be around 
3xl0 17 Mx. However, the inferred value is very probably af- 
fected by the limited spatial resolution of the TIP instrument. 
Therefore, the calculated total magnetic flux likely represents 
only a lower limit. 

The third row of Fig. [4] displays the temporal evolution of 
the magnetic field inclination which indicates significant differ- 
ences between the preexisting ambient field and the emerging 
one. While the former is vertical and stays almost constant in 
time, the latter shows a much more complex evolution. In its 
early stages, the inclination angles are close to 90°, correspond- 



ing to a horizontal magnetic field (see also the bottom panel of 
Fig. [5]l. Later on, the inclination on the edges of the feature be- 
comes more vertical. At 10:43:53 UT, the emerging feature has 
a clear loop-like structure with two opposite, almost vertical po- 
larities connected by the horizontal field. The trend for a de- 
crease of the inclination continues even after the linear polariza- 
tion signal has vanished; the remaining tracers of the loop foot- 
points become still more vertical and also separate further (see 
also the lower panel of Fig-EJ. In the final stage, the upper leg 
meets with a preexisting magnetic structure of opposite polarity. 
Reconnection of these two magnetic elements may thus explain 
the fast disappearance of the polarization signal at this location. 

The bottom row of Fig.|4]shows the magnetic azimuth angle. 
The emerging magnetic structure has an almost constant azimuth 
of ~75° during the whole time of its occurrence that is consistent 
with a rigid magnetic dipole. This again supports the proposed 
model of a magnetic loop. Due to the 180° azimuth ambiguity, 
we cannot resolve without additional information whether the 
loop has an n-like shape and is emerging, or an U-like shape and 
is submerging (CE07). 

To resolve the mentioned ambiguity between an emerging 
or submerging structure, we analyzed the Doppler shifts of the 
magnetic inversion component. The most relevant maps, show- 
ing the phase of the appearance of the structure under study, are 
displayed in Fig. [6] We found that the middle part of the appear- 
ing structure was blueshifted all the time of its occurrence; the 
velocity reached up to -2 km s _1 . This suggests an emergence of 
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Fig. 4. Temporal evolution during the flux emergence. The G-band images are shown in the 1 st row. Areas with enhanced net linear 
polarization are marked by green contours. The blue and red contours represent positive and negative circular polarization, respec- 
tively. The outer contours enclose regions where linear and circular polarization are greater than 0.13% and 0.11%, respectively, 
and the innermost contours represent a 0.16% level for both kinds of polarization. The other rows show the magnetic flux density 
(2 nd row), inclination (3 rd row), and azimuth (4 th row). The gray contours mark the same polarization levels as before. The times 
given on top correspond to moments when TIP passed across the occurring feature during the scanning. 



a magnetic loop from the lower atmospheric layers. Only later, 
when the footpoints of the loop become clearly visible, we could 
detect redshifts at the locations of the footpoints, corresponding 
to downflows of up to +2 km s _1 . This could signify that the ma- 
terial in the loop is draining down along the field lines, which 
is in agreement with the generic loop emergence scenarios. The 
surrounding magnetic structures are on the other hand predomi- 
nantly redshifted during the whole displayed time interval. 

Based on the observational findings discussed above, we 
conclude that the temporal evolution, the appearance of the 
linear and circular polarization signal, and the corresponding 
parameters of the magnetic field together with the 'magnetic' 
Doppler shifts derived by the inversion all suggest that the ob- 
served flux emergence takes place in the shape of a small loop 
that rises through the photosphere. 



5. Emerging loop in semiempirical simulation 

5. 1 . Geometrical model 

To investigate whether the observed event can be explained by 
flux emergence in the shape of a rising loop, we modeled an f]- 
shaped loop with the basic geometric parameters as in the obser- 
vations, and synthesized spectra of the IR lines for it. The loop 
was placed in a 3D box of around 3000kmxl000kmx500km, 



with a spatial sampling of 29 km in x and y, and 1 km in the z 
direction. 

The loop parameters derived from the observations were a 
horizontal separation speed of its footpoints of 2kms~ 1 , a hori- 
zontal extent of around 600 km at the beginning, and a life time 
of around 600 s, sampled on 6 maps of the repeated scanning. We 
then created a loop connecting two footpoints with a distance of 
L(k) =600 km + k x 180km (k = 1,2, 6) by prescribing a 2 nd 
order polynomial of the shape 



z(x) = h 



h 



apex 



apex 2 
X . 



(axL/2) 2 ' 



(1) 



where a = 0.6 and h apex is the location of the apex of the loop. 
The value of a influences the shape of the loop by stretching 
or compressing it. It was also used to take into account that 
the height scale considered for the later conversion to optical 
depth extended beyond km to z — -86.2 km, whereas the cross- 
section seen in the spectra rather comes from z > km. The loop 
parameter h apex was determined by trial-and-error to match the 
temporal evolution in the observations to h apex (k) = -65 km + 
fcxlOOkm, corresponding to a rise speed of around 1 kms~'. 

Figure|7] shows the resulting shape for the central axis of the 
loop on the six steps on the top. The inclination of the magnetic 
field to a vertical LOS (2 nd panel) was taken from the derivative 
of the location of the central axis. We determined an appropriate 
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Fig. 5. Cuts along the emerging feature showing the temporal 
evolution of its magnetic flux density (top) and inclination (bot- 
tom). 
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Fig. 6. 2D maps showing temporal evolution of the 'magnetic' 
Doppler shifts. The displayed values are cut at +2 km s . 



LOS velocity of the material inside the loop from two different 
contributions, the rise of the loop itself (blueshift in the middle) 
and a downflow along the field lines (redshift at the footpoints). 

The velocity component induced by the rise of the loop was 
determined from the vertical motion of the loop as defined by 
Eq. ([]]). We parametrized the loop for every step k by equidistant 
pieces U of length L = 1km along its axis (i = Q,..n max , where 
rimax was set by where the loop reached z = -1000 km), starting 
at the apex at each step, to obtain the positions x(k, If), z(k, If) in 
the x-z plane (see Fig. [8]). The velocity component for a vertical 
LOS at the position x(k + 1 , /,•) on the abscissa was then deter- 
mined as v z (x(k+ I, If)) = (z(k + 1, /;) - z(k, h))/dt, where dt cor- 
responds to the ~100 s cadence of the TIP maps. To obtain v z for 
the evenly spaced grid points used in the spectral synthesis, the 
velocity values v z (x(k + 1, /,•)) were then interpolated to equidis- 
tant 29km-spacing in x. The velocity ranges from -1 kms -1 at 
the apex to slightly higher velocities of about -1.4km s _1 near 
the footpoints due to the additional lateral motion (see 3 panel 
of Fig. |7]l. For the downflow along the field lines, we assumed 
a steady flow in free-fall conditions along the loop geometry 
of step k-b, similar to Mein et al. (1 19961 1 or Georgakilas et al. 
( 1990 their Eq. (11)). The free-fall condition, however, produces 
supersonic flows of v »7kms~ 1 near the footpoints. We thus 
have set a limit for the maximum allowed absolute velocity of 
4kms~', which leads to around +2kms _1 vertical velocity near 
the footpoints (bottom panel of Fig.|7]l. 



5.2. Magnetic and thermodynamical properties 

We modeled the loop as a circular flux tube around the central 
axis by using a field strength B with 



B(r) = 500Gxexp 



2cr 2 



(2) 



where r measures the distance from the central axis and cr was 
set to 100 km. 

This yields a flux tube with around 240 km FWHM and an 
effective radius of 300 km, at which distance the field strength 
drops to virtually zero. Integration of Eq. (f2|i over the radius 
yields a total magnetic flux of around 2.7xl0 17 Mx, which is 
slightly below the magnetic flux value derived from the inver- 
sion results. Figure|9] displays the final geometry by showing 
isosurfaces of 100G field strength for the six steps of the flux 
emergence. At the first step, most of the loop is below the for- 
mation height of the spectral lines considered, and thus induces 
only a small polarization signal in the spectra. 

We used the tabulated values of the HSRA model (Gingerich 
et al. 1 197 It to define the thermodynamic properties like temper- 
ature, gas density, or electron pressure needed for the spectral 
synthesis. We considered an optical depth range from logr = 
+ 1.4 to around -3, which covers the formation heights of the 
IR lines at 1.56 micron (Cabrera Solana et al. 2005). We did not 
consider the effect of the presence of the loop on the stratifica- 
tions (or vice versa) by enforcing hydrostatic equilibrium. 

To determine the magnetic field stratifications for a verti- 
cal LOS from the 3D topology, we used cuts through the 3D 
box along the z-axis. The field strength B(x,y,z) can then be 
derived from the distance of the point (x, y, z) from the closest 
piece of the central axis using Eq. ©. The stratifications of B 
can differ strongly from a Gaussian shape, depending on where 
along the loop the cut is located; near the footpoints, one obtains 
a Gaussian shape stretched in z, whereas close to the apex the 
Gaussian shape is preserved (Fig.fTUIl. 



Gomory et al.: Magnetic loop emergence within a granule 



7 





< -1-2 
E 

# _, - 3 t~ 

-1.4 h 

-i-5i_ 







[km] 



Fig. 7. Temporal development of the simulated loop in six time 
steps. From top to bottom: Central loop axis, magnetic field in- 
clination, LOS velocity component v z of the loop rise, free fall, 
and their addition to the final v,. 



5.3. Spectral synthesis and degradation of spatial resolution 

Using the derived stratifications for the thermodynamic variables 
and the magnetic field strength, we then synthesized spectra of 
the three spectral lines Fei 1564.8 nm, a weak blend of Fei at 
1564.7 nm, and Fe 1 1565.2 nm with the SIR code. The simulated 
spectra differ strongly from the observations in spatial resolu- 
tion. To improve the match to the observations, we degraded the 
spatial resolution by a convolution with a 2D Gaussian kernel of 
(X'56 FWHM, roughly corresponding to the diffraction limit of 
the VTT at IR wavelengths. We then calculated 2D maps of total 
linear and circular polarization L and \V\, respectively, and the 
polarity of the Stokes V signal analogously to the observations. 

FigureQT]shows these maps for both the original and the spa- 
tially degraded simulation, and the corresponding section from 
the observation, for a visualization of the effects of spatial reso- 
lution and spatial sampling. The display ranges of each quantity 
in Fig.[TT]are given in TableQ] For the observations, the value of 
the wavelength-integrated quantities cannot reach zero because 
of the noise. With the assumption that the lowermost value in 
the observations corresponds to a pure noise signal, the dynamic 
range of total linear polarization L in the observations matches 
fairly well that of the convolved simulation, whereas for total 
circular polarization |V| it is significantly smaller. The spatial 
smearing is sufficient to reduce the circular polarization signal 
in the first step to values below the threshold used for the ob- 
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Fig. 8. Calculation of velocity components (y x , v z ) from the lo- 
cation of the central axis in the steps k and k + 1 . 



servations (compare 7 s ', 7 th , and 13 th column). The elongated 
loop structure seen in L of the full-resolution spectra (top left) 
is lost after the spatial degradation; the strong linear polarization 
signal appears there in the shape of an only slightly elongated 
patch. The two polarities in the degraded synthetic polarity map 
are seen to touch each other persistently up to step 3, similar to 
the observed behavior. The effect of applying the Gaussian ker- 
nel is to first order already equivalent to sampling the degraded 
data with the pixel size of the observations (white rectangle in 
the bottom right panel, 0'.'45x0'.'17 pixel). We thus skipped an 
additional re-sampling with the TIP pixel size for the simulation 
spectra shown in the following. 

5.4. Direct comparison of observed and synthetic spectra 

For a direct comparison of observed and synthetic spectra, we 
took cuts along the central axis of the feature in the simulation. 
For the observations, we chose the slit spectra of a scan step 
near the beginning of each repeated map, since the structure is 
roughly oriented along the slit. FigureQ~2]shows the IQUV spec- 
tra for these six cuts. The observed spectra (at bottom) show 
some features that are naturally missing in the simulation (at 
top): the intensity and velocity pattern of the granulation, and 
additional polarization signal at the upper end of the cut corre- 
sponding to the preexisting patch of network fields (y ~ A"). 
For a comparison of the polarization level in the simulation and 
observation, we added an rms noise of 1 x 10 -3 to the simula- 
tion spectra of Stokes Q of the 1565.2 nm line and Stokes U of 
the 1564.8 nm line in the first three steps (white rectangle at top 
middle). The noise level is sufficient to suppress the polarization 
signal of the less Zeeman-sensitive line at 1565.2 nm as in the 
observed spectra. The maximal linear polarization amplitudes of 
the simulation and observation also roughly match (up to 1.5% 
of I c for Q and U). The Stokes V amplitude of the simulation 
exceeds the observations for the steps 5 and 6, where in the ob- 
servations the cancellation with the preexisting network patch 
may already have happened. The observed trend of polarization 
signal with time is reproduced in the simulation: a monotonical 
increase of V amplitude until step 5, increase and decrease for 
Q and U with maximal amplitude at step 4. The observed spec- 
tra still show a faint Q and U signal in step 6, only seen in the 
1564.8 nm line. Near the footpoints, the observed Stokes V sig- 
nal has a drop-like shape, with a decrease of the splitting and a 
displacement to the blue when moving towards the opposite po- 
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Fig. 9. Magnetic field topology during the flux emergence simu- 
lation. 

Table 1. Display ranges for Fig.fTTI 
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larity that also shows up in the simulation spectra. The observed 
Stokes Q signal in step 3 has the same curved shape in y as the 
synthetic spectra of steps 3 and 4, indicating the blueshift in the 
middle and the redshift near the footpoints. 

At step 4 (10:42:10UT), the loop-like structure is seen most 
clearly in the observations (cf. Figs. l4lor[TTt. We thus selected 
the profiles at the center of the loop and both footpoints of this 
step from observations and the simulation (short horizontal bars 
in Fig.n~2b for the direct comparison shown in Fig. [13] Besides 
the amplitude of Stokes V, the continuum intensity level and 
the Doppler shift in the lower footpoint, the match is surpris- 
ingly good, despite the simplifications used in generating the ge- 
ometric model. We thus conclude that the observed spectra are 
in good agreement with a rising magnetic flux tube of around 
250km diameter, 3xl0 I7 Mx magnetic flux, and a rise speed of 
1 kms -1 . 



6. Discussion and conclusions 

High-resolution surveys of large areas of quiet-Sun internet- 
work (Lites et al.|2008) and plage regions (Ishikawa et al. 2008 ) 
revealed widespread occurrences of relatively strong horizon- 
tal fields over granules, accompanied by vertical bipolar fields 
cospatial with up- and downflows (Martinez Gonzalez et al. 
2007). On the other hand, individual granules are often collo- 



cated with strong unipolar fields decoupled from granular out- 
flows (OR08, Bello Gonzalez et al. |2008l . 

CE07 reported for the first time the emergence of a magnetic 
loop inside of a granule whose footpoints drifted towards inter- 
granular lanes. Here we present new observational evidence for a 
similar event, whose onset manifests itself through enhanced net 
linear polarization within an evolved granule. Soon after, a pair 
of opposite circularly polarized patches appear symmetrically on 
each side of the linear polarization forming a small magnetic 
dipole bridged by a horizontal magnetic field. As long as both 
patches of circular polarization are visible, they are about of the 
same size and show a balanced magnetic flux of opposite sign 
(+2.9 and -3.1 x 10 17 Mx in step 4). This supports the conclu- 
sion made by Lamb et al. (2008 ) that the model of an asymmetric 
dipole emergence is very unlikely. 

The fast disappearance of the linear polarization signal after 
10:43:53 UT (step 5, see Figs.[TJ[4] orfTTT) is in striking contrast 
with its gradual growth seen from 10:38:47 UT to 10:43:53 UT. 
Can we understand it just by an upward extending loop reach- 
ing up to a higher atmospheric layer undetectable by the Fei 
lines as suggested in CE07? If we consider the loop top as a 
magnetic perturbation propagating upwards through otherwise 
non-magnetic layers, then a smooth decay of contribution and 
response functions with height would imply slower ceasing of 
the net linear polarization, perhaps detectable on the next frame. 
Since the contribution and response functions of Stokes Q and U 
to magnetic perturbations can be very complex (del Toro Iniesta 
2003), their detailed computations would help to understand the 
disappearance of linear polarization after 10:43:53 UT. But after 
10:43:53 UT the circular polarization signal of the footpoint that 
meets an opposite-polarity magnetic field also weakens rapidly. 
This could indicate magnetic reconnection between the emer- 
gent flux and the preexisting magnetic fields. If reconnection 
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Fig. 10. Magnetic field strength in side view in the x - z plane 
(top), and cut in the y — z plane across one of the footpoints (bot- 
tom). The location of the cut is marked by a dash-dotted vertical 
line in the upper plot. 
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Fig. 11. Top to bottom: Total linear polarization, circular po- 
larization, and polarity. Left: Full resolution of the simulation. 
Middle: After degradation. Right: Corresponding section of the 
FoV in the observations. The white rectangle in the bottom right 
panel indicates the size of a pixel in the observations. Contour 
lines trace strong linear polarization signal. Each panel and 
quantity is displayed with an individual range for better visibility 
(see text for details). Tickmarks are in arcsec. 









II 
I) 


m 


1 

1 


6 








l{ 




1 

I 


5 








II 






4 








u 








3 

l 

| 











+ 0.1% 


noise 




2 








Q 


+ 1.5% 


U 


±1.5% 


V 


1 

±5% 












6 








1, 




1 


5 








II 




• 


4 








(| 


1 


3 

1 












2 


1 j 




Q 


±1.5% 


u 


±1% 


V 


1 

±1.8% 


0.0 


0.2 0.4 


0.6 


0.0 0.2 


0.4 0.6 


0.0 0.2 


0.4 0.6 


0.0 0.2 


0.4 0.6 



* -1554.8 nrr 



Fig. 12. Comparison of observed {bottom) and simulated spectra 
(top). Left to right: Stokes IQUV. Bottom to top: Time steps 1 
to 6. The display ranges for QUV are given at bottom right in 
each panel. The vertical dashed line in Stokes / denotes the rest 
wavelength of the 1564.8 nm line. The short horizontal bars at 
the left in IQUV of step 4 denote the location of the profiles 
shown in Fig. [13] 
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thus is such a common feature, it could explain why OR08 found 
no trace of the flux emergence in the chromospheric layers. 

The disappearance of the polarization signal could, however, 
also be related to other causes. In the rising loop simulation, 
the linear polarization amplitude at the last step is again com- 
parable to the first one. Taking into account the noise present 
in the observations, the polarization signal could also be sim- 
ply reduced below the detection limit. Another possibility could 
be a slight spatial drift of the feature with time that would re- 
move it partly out of the observed FoV, because it already was 
located only in the first three steps of each repeated scan (see 
Fig-EJ- Moreover, the possible expansion of the magnetic loop 
with height (and time) could also affect the measurement of the 
linear signal above the detection threshold. From the comparison 
with the simulated loop, we conclude, however, that at least the 
ascent until maximal linear polarization signal is reached is fully 
traced by the observations. 

Our flux density estimate of ~450 G found in the emerging 
loop is comparable to the values reported in Martinez Gonzalez 
et al. (120071 1 and Ishikawa et al. (2008) for small-scale mag- 
netic loops in the solar internetwork and plage region, respec- 
tively; it significantly exceeds, however, the 10-100G given by 
Martinez Gonzalez & Bellot Rubio d20091 l. This could be due to 
the assumption of a magnetic filling factor of unity by the lat- 
ter authors. The loop footpoints separate with time by around 
2kms~', drifting in opposite directions towards intergranular 
lanes. Moreover, our estimate of the magnetic flux density shows 
that the field is stronger than the typical equipartition value of 
400 G (Cheung et al. 120071 Ishikawa et al. 2008). This suggests 
that both granular outflow and relaxation of the magnetic tension 
of the loop may drive the footpoints. 

How does our loop compare to the one reported in CE07? 
They occupy similar locations and enclose too little magnetic 
flux to leave any trace on the underlying granules or ambient G- 
band bright points. On the other hand, our loop displays a higher 
magnetic flux density and a longer lifetime. 
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Another difference to the loop reported in CE07 is that in our 
case one loop footpoint meets the ambient field of opposite po- 
larity, which is followed by the fast disappearance of the loop, 
in contrast to its more gradual emergence. This is very symp- 
tomatic of field cancellation by magnetic reconnection, as re- 
ported in Cheung et al (120081 1 and Guglielmino et al. ( 120081 . The 
magnetic reconnection likely drags loop remnants very rapidly 
out of the line-forming domain less than 12 min after the first 
trace of the loop top was detected. 

The magnetic loop emergence presented here is probably just 
a particular event in a multitude of similar ones occurring con- 
tinuously in the solar internetwork. Simulations of small-scale 
magnetoconvection by, e.g., Stein & Nordlund (2006) show ex- 
amples of granular-sized magnetic loops with similar behavior. 
They probably build up a sub-canopy of horizontal fields over 
granules (Wedemeyer-Bohm et al. 2009). Surface dynamo sim- 
ulations by Schussler & Vogler (2008) and Steiner et al. (2008 ) 
clearly show that this new component of the internetwork is a 
direct consequence of magnetic flux expulsion from the granu- 
lar interior, not only to the intergranular lanes, but also to the 
upper photosphere. Cheung et al. (2007 2008 ) found in numeri- 
cal simulations that emerging flux regions also contribute to the 
flux accumulation and its maintenance above granular upflows 
in the magnetic inversion layer. Our results show one particular 
example of this local process albeit terminated, probably due to 
magnetic reconnection. Better temporal and spatial resolution is 
needed for more detailed information on the process than our 
observations provide. 
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